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In a previous paper, we computed expressions for the Detweiler- Whiting singular field of point 
scalar, electromagnetic and gravitational charges following a geodesic of the Schwarzschild spacetime. 
We now extend this to the case of equatorial orbits in Kerr spacetime, using coordinate and covariant 
approaches to compute expansions of the singular field in scalar, electromagnetic and gravitational 
cases. As an application, we give the calculation of previously unknown mode-sum regularisation 
parameters. We also propose a new application of high order approximations to the singular field, 
showing how they may be used to compute m-mode regularization parameters for use in the m-mode 
effective source approach to self-force calculations. 



I. INTRODUCTION 

The two-body problem in general relativity is a longstanding open problem going back to work by Einstein himself. 
With recent advances in gravitational wave detector technology, this age-old problem has been given a new lease of 
life. Some of the key sources expected to be seen by both space and ground based gravitational wave detectors are 
black hole binaries. Accurate models of black hole binaries are required for their successful detection by gravitational 
wave detectors. This development is today motivating numerical, analytical and experimental relativists to work 
together with the goal of producing models of the inspiral and merger of binary black hole systems. 

In modelling black hole binaries, it is widely accepted that for the scenario of an extreme mass ratio inspiral (EMRI), 
the self-force approach is the model of choice. EMRIs are expected to be seen by space-based detectors such as 
NGO/eLISA pQ. Although NGO/eLISA has been recently postponed, the gravitational wave research community are 
confident in its inevitable flight. In the meantime, recent research has shown the applicability of self-force calculations 
to other binary black hole configurations [H [3] , extending the application of self-force to ground-based detectors such 
as LIGO and VIRGO. 

Within the self-force approach, one perturbatively solves for the motion of a small body in the background of a 
massive black hole. Formal derivations of the equations of motion of a small body moving in a curved spacetime have 
settled on the idea of a well-defined singular-regular split of the retarded field generated by the body [3HH] ■ Several 
practical self-force computation strategies have developed from these formal derivations, all of which are based on 
the now-justified assumption that the use of a distributional source is acceptable at first perturbative order. These 
strategies broadly fall into three categories: the mode-sum approach p~3j [14], the effective source approach \TE\ [16] 
and Green function approaches [TTJ, HH] ■ The key to all three approaches is the subtraction of an appropriate singular 
component from the retarded field to leave a finite regular field which is solely responsible for the self-force. This 
singular component must have the same singular structure as the full retarded field in the vicinity of the body and 
must not contribute to the self-force (or its contribution must be well known and can be corrected for). There are 
many choices for a singular field that satisfies these criteria, although not all choices are equal. Detweiler and Whiting 
[T5] identified a particularly appropriate choice. Through a Green function decomposition, they defined a singular 
field which not only satisfies the above two criteria, but also has the property that when it is subtracted from the full 
retarded field, it leaves a regularized field which is a solution to the homogeneous wave equation. Extensions of this 
idea of a singular-regular split to extended charge distributions [5J [9] , to second perturbative order ptJHM] and fully 
non-perturbative contexts [25] have recently been developed. 

In a previous paper [26] (from now on referred to as Paper I), we focussed our calculations on the Schwarzschild 
spacetime representing a non-rotating black hole. Although this is a possible physical scenario, it is believed that a 
more astrophysically realistic or probable situation would be that of a Kerr or rotating black hole spacetime. One 
of the primary goals of the self-force community is, therefore, the successful calculation of the self-force in Kerr 
spacetime, with particular emphasis on the gravitational case. To this end, we now adapt our previous work [26j to 
the Kerr spacetime. 
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In Paper I, we developed approaches to computing highly accurate approximations to the Detweiler- Whiting singular 
field of point scalar and electromagnetic charges as well as that of a point mass. This was achieved through high-order 
series expansions in a parameter e, which acts as a measure of distance from the particle's world-line. We also derived 
explicit expressions for the case of geodesic motion in Schwarzschild spacetime. In this paper, we extend this analysis 
to the case of eccentric, equatorial orbits in the Kerr spacetime. We find that all of the methods developed in Paper I 
may be applied to the Kerr case with little modification. Nevertheless, the detailed expressions arc significantly more 
complicated in the Kerr case. Since our method is largely the same as that used in Paper I, we direct the reader there 
for full details and give here only the expressions which differ. 

As applications of our high-order expansions of the singular field, in Paper I we derived expression which may 
be used to improve the accuracy of both the mode-sum and effective source approaches. Since the effective source 
approach requires that the source be evaluated in an extended region around the world-line, numerical evaluation can 
be time consuming, in particular, when using high order expansions such as the ones produced in this paper and in 
Paper I. Existing calculations have settled on expansions of the singular field to 0(e 2 ) to be a particular 'sweet spot' 
for these calculations - up to this order the increase in complexity of the singular field and corresponding effective 
source is rewarded with an increase in accuracy. However, expansions above this order may well slow the calculations 
down to such a degree that the extra orders offer more of an hindrance than a help. In this paper we propose one 
solution to this problem, allowing most of the benefit to be reaped from high-order expansions without the need 
for using increasingly complicated high-order expansions in numerical evolutions. This idea makes use of the m- 
mode scheme, developed by Barack and Golbourn |15j . which decomposes the retarded field and effective source into 
azimuthal modes; the resulting conservation of axial symmetry makes the scheme well suited to the Kerr space-time. 
By carrying out m-mode effective source calculations with an effective source accurate to some order, say 0(e 2 ), one 
can obtain numerical values for the m-modes of the self-force which converges polynomially with 1/m. Our technique 
then makes use of our higher terms of the singular field (those above 0(e 2 )) to obtain a faster convergence of this 
m-mode sum and hence assist in the production of highly accurate values of the self-force. 

The layout of this paper is as follows. In Sec. [TTJ we use coordinate expansions to derive high-order regularization 
parameters for use in the mode-sum method. In doing so, we give the next two previously unknown non-zero 



regularization parameters in scalar, electromagnetic and gravitational cases. In Sec. Ill we propose a new application of 
high-order coordinate expansions of the singular field, showing how they may be used to derive m-mode regularization 
parameters for use in the m-mode effective source approach. In Sec. |IV| we summarize our results and discuss further 
prospects for their application. 

Throughout this paper, we use units in which G = c = 1 and adopt the sign conventions of [57]- We denote 
symmetrization of indices using parenthesis (e.g. (ab)), anti-symmetrization using square brackets (e.g. [ah]), and 
exclude indices from (anti-)symmctrization by surrounding them by vertical bars (e.g. (a|6|c), [a|6|c]). We denote 
pairwise (anti-)symmetrization using an overbar, e.g. R^—^ = \{R a bcd + Rcdab)- Capital letters are used to denote 
the spinorial/tensorial indices appropriate to the field being considered. In many of our calculations, we have several 
spacetime points to be considered. Our convention is that: 

• the point x refers to the point where the field is evaluated 

• the point x refers to an arbitrary point on the world-line 

In computing expansions, we use e as an expansion parameter to denote the fundamental scale of separation, so that 
x — x s» 0(e). Where tensors are to be evaluated at these points, we decorate their indices appropriately using ", e.g. 
T a and T a refer to tensors at x and x, respectively. 



II. /-MODE REGULARIZATION 



One of the most successful self-force computation approaches to the date is the mode-sum scheme of Barack and 
Ori [T31[T1]; the majority of existing calculations are based on it in one form or another [28145 1| . The basic idea is to 
decompose the retarded field into spherical harmonic modes which are continuous and finite in general for the scalar 
case and in Lorenz gauge for the electromagnetic and gravitational cases. The spherical symmetry of the Schwarzschild 
spacetime makes this decomposition into spherical harmonic modes a natural choice. In Kerr spacetime, despite there 
being more natural choices (such as a decomposition into spheroidal harmonics), a decomposition of the singular field 
into spherical harmonic modes has been shown to be of practical use in computing the scalar self- force [44 , 52 . While 
similar approaches have yet to be attempted in electromagnetic or gravitational cases, it seems likely that they are 
at least possible in principle. 

A key component of the mode-sum calculation involves the subtraction of the so-called regularization parameters - 
analytically derived expressions which render the formally divergent sum over spherical harmonic modes finite. In this 
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section, we derive these parameters from our singular field expressions and show how they may be used to compute 
the self-force with unprecedented accuracy. 



A. Mode Sum Concept 

The self-force for scaler, electromagnetic and gravitational cases, can be written generically as 

F a =p a Alfi f K) (2.1) 

where 

V(k) = yfrct) - P(S) ( 2 - 2 ) 

is the regularised field and p a a{x) is a tensor at x which depends on the type of charge. We can therefore rewrite the 
self-force as 

F a =p a Ai pf let) -p a Ai pf s) . (2.3) 
Carrying out a spherical harmonic decomposition on the field, 

CO I 

^et)/(S)=E E PMret)/(S)*W0,0), ( 2 - 4 ) 
i=o m=-l 

allows the self-force to be rewritten as, 

oo / 



pa = E E (p a ^L(rct) ~ P a A<pL {S) ) Y lm (9 , M- (2.5) 

1=0 m =-l 

Defining the I component of the retarded or singular self-force to be, 

l 

^T(ret)/(S) = P a A ^ ^(ret)/(S)^m(^0,^o) ) (2-6) 
m— — I 

the self-force can be expressed as 

oo 

^E^Vt)-^))- (2.7) 



(=0 



It is the last term on the right that we calculate in this section for each of scalar, electromagnetic and gravitational 
cases in Kerr spacetime. 

Our explicit expression for the Z-modes of the singular self-force in Kerr spacetime is written as an expansion about 
the world- line point x, that is 

^"(S) = F l aH] (r 0) to) e- 2 + F l a[0] (r , to) ^ + F l a[2] (r , to) e 1 

+ F l a[i] (r , h) e 3 + F l a[6] (r , t ) e 5 + . . . , (2.8) 

where we are missing even powers of e above —2, as these are zero - this will be shown to be the case later in this 
section. When summed over I, the contribution of FjjL (r ,i ), and higher terms, to the self-force is zero. However, if 

we ignore these higher terms in the approximation of ¥?^„(s)> then the approximation for f^m) is only C 1 , causing the 
sum over I to be polynomially, rather than exponentially convergent in l/l. Therefore, despite these terms having zero 
contribution to the self-force, when it comes to numerically calculating the self-force using a finite number of I modes, 
the inclusion of the higher order terms dramatically reduces the number of modes required and hence computation 
time. For this reason, every extra term or regularisation parameter that can be calculated is important. 
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B. Rotated Coordinates 



In order to obtain expressions which are readily written as mode-sums, previous calculations [131 1301 132] found it 
useful to work in a rotated coordinate frame. In Paper I, we found it most efficient to carry out this rotation prior 
to doing any calculations; this also holds in the Kerr case. To this end, we introduce coordinates on the 2-sphere at 
x in the form 



W\ = 2 sin J cos /3, w-i = 2 sin ^— ^ sin/3, 
where a and /3 are rotated angular coordinates given by 



(2.9) 



sin cos <j> = cos a, (2-10) 
sin 8 sin <fi = sin a cos /3, (2-11) 
cos9 = sin/3sina. (2-12) 

The Kerr metric in these coordinates (with x chosen to lie on the equator, ie. a = 0) is given by the line element 



ds 2 



8Mr 



Ar 2 + a 2 w 2 (4 — w\ — wj) 



- 1 



dt 2 



4r 2 + c?w\ (4 — w\ — uj|) 



4 (r 2 - 2Mr + a 2 ) 



2dtdw\ 



-2aMr [8 - w\ (6 - w\ - wj)] 



•)\ — w\ [Ar 2 + a 2 w 2 (4 — w\ — w?,)] 
1 



2dtdw? 



dr 2 

2aMrw\W2 (6 — w 2 — w|) 



74" 



w\ [Ar 2 + a 2 w 2 (4 — w\ — w. 



A(4-w 2 -w 2 )[A-w 2 (4- 



\g Wl w 1 dw 1 + 2g WlW2 dwidw 2 + . 



(2.13) 



where 



w l w 2 + 



,2„,,2 



a )] + [8 - w\ (6 - w\ 
xw 2 (w\ + 2w\ - 4) [Ar 2 + a 2 wl (A-w 2 - w%)] 



2 i 2 

r + a 



2Ma 2 r 



w 2 (A- 
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Ar 2 + a 2 w\ (4 — w\ — w 2 ) 



+ W1W2 (6 — w\ — iuf ) [8 — w 2 (6 



IV. 



D] 



a 2 + 2Mo 2 r 



(4 



4r 2 + a 2 w\ (A — w 2 — w 2 ) 



gw 2 w 2 = {A-w\- wlf [Ar 2 + a 2 wl{A - w\ - wfj] 



1 1 

w l w 2 



(6 



IT 



1 2 + 2Ma 2 r 



(4 



^1 



4r 2 + a 2 w 2 (A — w\ — w 2 ) 



(2.14) 



As in the Schwarzschild case, this algebraic form has an advantage over its trigonometric counterpart in computer 
algebra programmes where trigonometric functions tend to slow calculations down. Despite the apparent complexity 
of the Kerr metric in this form, calculations of the regularization parameters are more efficient than using Boyer- 
Lindquist coordinates and rotating the resulting complicated expressions. 



C. Mode Decomposition 

Having calculated the singular field using the Kerr metric in the above form and using the methods described in 
Paper I, it is straightforward to calculate the singular component of the self- force, F a , for scalar, electromagnetic 
and gravitational cases using Eq. (2.1 ) with the singular field substituted for the regular field. We obtain a multipole 
decomposition of F a by writing 

00 1 

F a (r, t, a,P) = J2 F im (r, *) Y lm (a, 0) , (2.15) 

1=0 m=-l 

where Y lm (9, <fi) are scalar spherical harmonics, and accordingly 



Ff m (r, t)= F a (r, t, a, (3) Y? m (a, /?) dSl. (2.16) 
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To calculate the /-mode contribution at x = (to, ro, ao, /3o)j we have 



Ft (r , h) = lim V F t a m (r + Ar, to) (a , A>) ■ 

Ar— >0 ' — ' 



(2.17) 



With the particle on the pole in the rotated coordinate system, Yj m (eta — 0, fto) — for all ra/0. This also allows 
us, without loss of generality, to take fto = 0. Taking a 0j fto an d m all to be equal to zero in Eq. (2.17) gives 



if (r , t ) = lim \ 2l -^F l l (r Q + Ar, t ) 



Ar^O V 47T 
2/ + 1 

47T A 



lim / F a (r + Ar,t , a, /3) Pi (cos a) dn. 



(2.18) 



Using the methods of Paper I, a coordinate expansion of the singular self-force F a (r, t, a, ft) may be written in the 
form 



°° Da(3n-2) 

F>,t,a,/3) = £-^r- £ «- 3 , 

n— 1 



(2.19) 



where 5 a ( fe ) = b a aia2 _ ak (x)Ax ai Ax a2 . . . Ax a <*. In using Eq. ( |2.19| to determine the regularization parameters, we 
only need to take the sum to the appropriate order: n = 1 for ~Fu_-n, n = 2 for i^/fop e ^ c - Explicitly, in our rotated 
coordinates 



p(r,t,a,PY 



Ar 2 r [r (a 2 E 2 - L 2 ) + 2M(L - aE) 2 + E 2 r 3 ] 



(a 2 -2Mr +r 2 Y 



At 



Awf 



2a 2 M 



a 2 +L 2 + rl 



2ArAwiLr 2 r 



Aw x ( 



At 2 E 



ro 



- 2EL) - 



2ArEr 2 r 
a 2 - 2Mr a + r 2 , 



2AI 
ro 



1 ] +Aw^ , 



(2.20) 



r "/ a 2 -2Mr + r 2 

where the a, ft dependence is contained exclusively in Aw\ and Aw2- Here, E — —Ut and L = are the energy 
per unit mass and angular momentum along the axis of symmetry, respectively. In particular, taking t = to (At = 0) 
allows us to write 

2 Ar 2 r [Erp (a 2 + rg) + 2aM(aE - L)} 2 

p(r,t ,a,p) = ^5 V AwnTr. 

(a 2 — 2Mr + rg) [r (a 2 + L 2 ) + 2a 2 M + r%] 



2a 2 M 2 r2 2 
+ a + L +r 2 



ArLrl 



o r o 



[a 2 - 2Mr + rg) (2a 2 M + a 2 r + L 2 r + rg) 



For the mode-sum decomposition, it is favourable to work with p Q (a, ft) 2 = p (r , to, a, ft) 2 in the form 

po {a, ft) 2 = 2(1- cos a) C 2 (l - k sin ft 2 ) 



This can be achieved by rewriting Eq. (2.21 1 with Ar — > as 



where 



and rearranging to give 



Pa (a,ft) 2 = C 2 Aw 2 1+ r 2 Aw 2 2 , 



9 o 9 2(2 J\d 9 

C 2 = L 2 + rl + + a 2 

ro 



po (a, ft) 2 = 2(1 — cos a) ( 2 



C 2 - r 2 



sin ft 2 



(2.21) 

(2.22) 
(2.23) 
(2.24) 

(2.25) 



which is equivalent to Eq. (2.22) with k = - ^ 2 r ° . Defining x(ft) — 1 — ksm 2 ft, we can rewrite our Aid's in the 
alternate forms 



Aw\ = 2 (1 - cos a) cos 2 ft = ^- cos 2 /3 

C x 

o 2 

Aw 2 = 2 (1 - cos a) sin 2 ft = ^- sin 2 /? 

C 2 x 



Po- 



IC 2 ~r 2 ) x 

2 



(fe-(i-x)) 



Po 

(C 2 -r 2 ) X 



(i-x). 



(2.26) 
(2.27) 



() 



Suppose, for the moment, that we may take the limit in Eq. ( 2.18 1 through the integral sign, then using our alternate 
forms we have 



lim 



j^a(3n— 2) 



n— 3 



b a 



(r )Aw il Aw i2 . ..Aw' 



jn,— 3 



Ar->0 p 



2n+l 



Pa 



2n+l 



(2.28) 



In [53], it was shown that the integral and limit in Eq. (2.18) are indeed interchangeable for all orders except the 
leading order, n = 1 term, where the limiting l//?o 3 would not be integrable. Thus we find the singular self-force now 
has the form 



Ft(r ,t ) 



2Z + 1 

47T 



B a ^ (r,t ,q,ft) 
Ar^oJ p 3 (r,t ,a,ft) 



2 lim 



oo - 

P, (cos a) dQ + ^ / Po n ~ 3 c a (n) (ro,x) Pi (cos a) 



zfyf.!] (r„, to) e" 2 + ^fo] (ro, *b) e" 1 + Ffa ( r o, to) e 1 + Ffa (r , to) e 3 + Ffa (r , to) e 5 + 



(2.29) 



where the ft dependence in the c? n j's are hidden in x, while the a and ft dependence of F a (r, to, a, ft) is hidden in both 
the p's and c? n )'s. Note here that we use the convention that a subscript in square brackets denotes the term which 
will contribute at that order in X/l. Furthermore, the integrand in the summation is odd or even under Awi — ¥ —Awi 
according to whether n (and so 3n — 2) is odd or even. As a result only the even terms are non- vanishing, while 
Ffa (r , t ) = Ffa (r , t ) = JO] (r 0) to) = 0, etc. 

Some care is required in order to obtain easily integrable expressions in the case of eccentric orbits. We use the 
approach of previous methods [30J 1321 H3J [M] (and also employed in Paper I) by redefining our Aw% coordinate in 
such a way that the cross-terms involving ArAwi in po vanish. That is, we make the replacement Awi — > Awi +cAr, 
where c is given by 



-Lr& 



This allows us to write 



(a 2 - 2Mr + r 2 ) (2a 2 M + a 2 r + L 2 r + rg) ' 



p (r, t , a, ft) 2 = v 2 Ar 2 + ( 2 Aw 2 + r 2 Aw 2 



(2.30) 



= v 2 Ar 2 + 2xC 2 (1 - cos a) 
where v is an expression involving ro, a, E and L. This can be easily rearranged to give 

p (r, t , a, ft)" 3 = C 3 (2*r 3/2 (5 2 + 1 - cos a)~ 3/2 

= C 3 (2 X )- 3/2 ^^- 3/2 (<5)^ (cos a) , 



;=o 



where 



S 2 = 



K 2 x 



and A x 5 (6) = 



21 + 1 



(2.31) 



(2.32) 



(2.33) 



Here A x 2 (S) is derived from the generating function of the Legendre polynomials as shown in Eq. (D12) of [3D]. We 
can now express p (r, to, a, ft) 3 as 



p(r, t , a, ft)" 



C 2 vxv Ar 2 l=Q 



^(Z + i)P(cosa) 



(2.34) 



Bringing this result into our expression for Fj^^ (ro,t ) from Eq. (2.29) and integrating over a gives 



(»*>*>) = ^(1 + 1)^ 



Ar^O c, 2 v\J 'Ar 2 

K Ar 1 
= Z+i lim — 

V 2; Ar->0 (^vy/Ej* 27T 

C^o 



B a ^Pi (cos a) 



X^dft 



(2.35) 
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where the first equality takes advantage of the orthogonal nature of the Pi (cos a) and the last equality comes from 
taking the limit as Ar — > and noting from Appendix C of [30 that the integral is a special case of the hypergeometric 
functions given by 

i/x^-ZftH*)-^-^. (2-36) 



B a ^ and also now carry a tilde to signify that they are not the exact same i? a W and b a a from Eq. (2.19); the 
tilde reflects the fact that they hav e also undergone the coordinate shift Aw\ ~ > Aw± + cAr. 

In the higher order terms in Eq. (2.29), we may immediately work with p 2 = 2x£ 2 (l — cos a) so, 

po (r ,t (h a,(3) n = C [2 X (1 - cosa)]" /2 

= C (2 X ) n/2 ^^ /2 (0)fl (cos a) , (2.37) 

1=0 

where A l 2 (0) = v2 from the generating function of the Legendre polynomials and, as given in Appendix D of 30J, 
for (n + l)/2 e N 

W2fn) = V n/2 (2/ + 1) 

; 1 j {2l-n){2l-n + 2)...{2l + n)(2l + n + 2y y ' ; 

(2.39) 

V n/2 ={-l) [n+1),2 2 1+n l 2 (nil) 2 . 



where 



(2.40) 

In this case the angular integrals involve 

1 f ^ U- n/2 (P))=^(lil,k) (2.41) 



= jr_x(fc)= 2 F 1 (-i ,l;l;k) =~£, (2.43) 



2tt J xW) n/2 \ I V 2 ' 2 ' 

where (n + l)/2 e N U {0}. The resulting equations can then be tidied up using the following special cases of 
hypergeometric functions 

X~ i )=J r i (k) =2^(^51;*) =1%, (2-42) 

1 

2' 2' 
where 

K= / (1 - fcsin 2 /3r 1/2 d/3, £ = / (1 - fc sin 2 f3) 1/2 d/3 (2.44) 
Jo Jo 

are complete elliptic integrals of the first and second kinds, respectively. All other powers of x can be integrated 
to give hypergeometric functions that can then be manipulated to be one of the above by the use of the recurrence 
relation in Eq. (15.2.10) of [55], 

= -^IS^-iik) + ^—^r^r^Mk)- (2-45) 
p(fc-l) P(K-l) 

In the next sections, we give the results of applying this calculation to each of scalar, electromagnetic and gravita- 
tional cases in turn. In doing so, we omit the explicit dependence on I which in each case is 

pa 

F l{ _ x] -(2l + l)F H] , F l[0] -F [0] , F l[2] - {2l _ 1){2l + 3 y 
l[4] (2i-3)(21-l)(2i + 3)(2i + 5)' 

pa 

pa _ [6] ley AC\ 

m (2Z - 5)(2Z - 3)(2Z - l)(2f + 3)(2Z + 5)(2Z + 7) ' 1 ' ; 
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D. Scalar i-raode Regularization Parameters 



In the scalar case, the singular part of the self-force is give by, 

F a =d a & s \ 

where <!> < > s ) is the scalar singular field. The scalar regularization parameters are then given by, 

r r sgnAr 



t[ " 1] ~ r (a 2 + L 2 ) + 2a 2 M + r 3 ' 
sgnAr [Er Q (a 2 + r 2 ) + 2aM(aE - L)) 



(2.47) 



(a 2 - 2Mr + r 2 ) (r (a 2 + L 2 ) + 2a 2 M + rg) ' 
F^.!] = 0, - 0, 



(2.48) 



) 



3/2 



(2.49) 



where 



F t £ [Q] = AaLM (Aa 4 M 2 + 2a 4 Mr + 2a 2 L 2 Mr - a 2 Mr 3 - a 2 r 4 - L 2 r 4 ) 

+ E(- 12a 6 M 3 - Wa 6 M 2 r - 7a 6 Mr 2 - a 6 r 3 - 4a 4 L 2 M 2 r - 6a 4 L 2 Mrl - 2a 4 L 2 r 3 - 6a 4 M 2 r 3 - 5a 4 Mr 4 
- a 4 r 5 + a 2 L 4 Mr 2 - a 2 L 4 r 3 - 5a 2 L 2 Mr 4 - 3a 2 L 2 r 5 - 2L 4 r 5 ), 
F^ 0] = -2aLM(2a 4 M 2 - a 4 Mr - a 4 r 2 - a 2 L 2 Mr - 2a 2 L 2 r 2 - 2a 2 Mr% - 2a 2 r\ - L 4 r 2 - 2L 2 r 4 ) 

+ E(4a 6 M 3 + 4a 6 M 2 r + a 6 Mr 2 - 2a 4 £ 2 M 2 r - a 4 L 2 Mr 2 + 2a 4 M 2 r 3 + a 4 Mr 4 - 2a 2 L 4 Mr 2 + a 2 L 2 Mr 4 
+ a 2 L 2 rl + L 4 rl), 



Ff [0] = ( - 24a 8 M 3 r - 32a 8 M 2 r 2 - 14a 8 Mr 3 - 2a 8 r 4 + 24a 6 L 2 M 4 + 12a 6 L 2 M 3 r - 30a 6 L 2 M 2 r 2 - 27a 6 L 2 Mr 3 

- 6a 6 L 2 r 4 + 48a 6 M 4 r 2 + 40a 6 M 3 r 3 - 16a 6 M 2 r 4 - 20a 6 Mr* - 4a 6 r 6 + 8a 4 L 4 M 3 r - \2a 4 L 4 Mr 3 

- 6a 4 L 4 r 4 + 36a 4 L 2 M 3 r^ + 12a 4 L 2 M 2 r 4 - 21a 4 L 2 Mr i - 9a 4 L 2 r% + 24a 4 M 3 ^ + 8a 4 M 2 ^ - 6a 4 Mr^ 

- 2a 4 r 8 - 2a 2 L 6 M 2 r 2 + a 2 L 6 Mr 3 - 2a 2 L 6 r 4 + 6a 2 L 4 M 2 r 4 + a 2 L 4 Mr 5 - 6a 2 L 4 r 6 + \2a 2 L 2 M 2 r% 

- 3a 2 L 2 r% + 2L^Mr% - i 6 r 6 + 2L 4 Mr 7 - L 4 r 8 ) 

- 2aELM(24a 6 M 3 + 28a 6 M 2 r + 10a 6 Mr 2 + a 6 r 3 + 8a 4 L 2 M 2 r + 8a 4 L 2 Mr 2 + 2a 4 L 2 r 3 - Aa 4 Mr 4 - 2a 4 r^ 

- 2a 2 i 4 Mr 2 + a 2 L 4 r 3 - 2a 2 L 2 Mr 4 - a 2 L 2 r 5 - 6a 2 Mr 6 - 3a 2 r 7 + L 4 r 5 ~ 3L 2 r 7 ) 

+ E 2 (2a 2 M + a 2 r + r 3 ) (l2a 6 M 3 + 16a 6 M 2 r + 7a 6 Mr 2 + a 6 r 3 + 4a 4 L 2 M 2 r + 6a 4 L 2 Mr 2 + 2a 4 L 2 r 3 
+ 6a 4 M 2 r 3 + 5a 4 Mr 4 + a 4 r 5 Q - a 2 L 4 Mr 2 + a 2 L 4 r 3 + 5a 2 L 2 Mr 4 + 3a 2 L 2 r 5 + 2L 4 r 5 ), 
Ff m = (8a 8 M 3 r + 8a 8 M 2 r 2 + 2a 8 Afr 3 - 8a 6 L 2 M 4 + 4a 6 L 2 M 3 r Q + 12a e L 2 M 2 r 2 + 4a 6 L 2 Mr 3 - 16a 6 M 4 r 2 

- 8a 6 M 3 r 3 + 8a 6 M 2 r 4 + 4a 6 Mrl + 4a 4 L 4 M 3 r + 8a 4 i 4 M 2 r 2 + 2a 4 i 4 Mr 3 + 8a 4 L 2 M 3 r 3 + 12a 4 L 2 M 2 r 4 
+ 2a 4 i 2 Mr 5 - a 4 L 2 r e - 8a 4 M 3 r 5 + 2a 4 Mr 7 Q + 4a 2 L 6 M 2 r 2 + 16a 2 L 4 M 2 r 4 - 2a 2 L 4 r 6 + 4a 2 L 2 M 2 r 6 

- a 2 L 2 r 8 + 2L?Mrl - L 6 r 6 + 2L 4 Mr 7 Q - i 4 r 8 ) 

+ 2aELM(8a 6 M 3 + 4a 6 M 2 r - 2a 6 Mr 2 - a 6 r 3 - 4a 4 L 2 M 2 r - 6a 4 L 2 Mr 2 - 2a 4 L 2 r 3 - 8a 4 M 2 r 3 

- 12a 4 Mr 4 - 4a 4 r^ - 4a 2 L 4 Mr 2 - a 2 L 4 r 3 - l0a 2 L 2 Mr 4 - ha 2 L 2 r% - 6a 2 Afr 6 - 3a 2 r 7 - L 4 r% - 3L 2 r 7 ) 

- E 2 (2a 2 M + a 2 r + rg) (4a 6 M 3 + 4a 6 M 2 r + a 6 Mr 2 - 2a 4 L 2 M 2 r - a 4 L 2 Mrl + 2a 4 M 2 r 3 l + a 4 Mr\ 

- 2a 2 L 4 Mr 2 + a 2 L 2 Mr 4 + a 2 L 2 r 5 + L 4 r 5 ), 



3 (2a 2 M + a 2 r Q + L 2 r Q ) 2 (r 2 + L 2 + 2 -^ + a 2 ) (r 2 - 2Mr + a 2 ) 




(2.50) 



where 
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F e[0] = 0, (2.51) 



F ck[0] - ~. "172 ( F r[0] £ + F m^) > (2.52) 



irr {2a 2 M + a 2 r Q + L 2 r f (r 2 + L 2 



2 j_ 2a' 2 M j_ fl2 



where 



F r £ [01 = -2a 4 A/ 2 - a 4 Afr - a 2 L 2 Mr - 4a 2 Mr 3 - a 2 r% - L 2 r%, 



■ r[0] 



The regularisation parameters F a [ 2 ] and -F a [4] are too large for paper format and have instead been made available 
electronically [56] . For the reader to get an understanding of the form and size of these expressions, we have included 
F r [ 2 ] for a circular orbit. This is given by, 



F r [2] 



2a^jMr% - ZMr% 



r ° \ (a^M + r 3/2 ) (a 2 + r (r - 2M)) 

(2aVM + V^(r - 3M)) ^ (f* 2] £ + Fj^K 



(cfiM + 2aV Mr o + fl2r o {~2M 2 + Mr + r 2 ) - 4aAf 3 / 2 r 5/2 + Mr^j 



3 ' 



(2.53) 



(2.54) 



where 



PC _ 



(r 3/2 + a^Af) (3M 3 r 33/2 - 36M 4 r 31/2 + 9a 2 Af 2 r 31/2 - 78aM 7 / 2 r 15 + 93M 5 r Q 29/2 + 27a 2 Af 3 r 29/2 

+ 9a 4 Afr 29/2 + 744aM 9 / 2 r 14 + 44a 3 A/ 5 / 2 r 14 + 3a 6 r 27/2 - 60Af 6 r 27/2 - 1068a 2 Af 4 r 27/2 - 120a 4 Af 2 r 27/2 

- 1890aAf 11/2 ^ 3 + 60a 3 Af 7/2 r^ 3 + 18a 5 M 3/2 r^ 3 + 3681a 2 AT 5 r 25/2 + 913a 4 M 3 r 25/2 + 15a 6 Mr 25/2 
+ 1332aAf 13/2 r 1 , 2 - 704a 3 Af 9/2 r^ 2 - 404a 5 M 5/2 r 1 , 2 + 24a 7 VMr 1 1 2 + 6a 8 r 23/2 - 4251a 2 A/ 6 r 23/2 

- 3846a 4 Af 4 r 23/2 - 357a 6 M 2 r 23/2 + 2472a 3 M 11/2 r^ + 3434a 5 Af 7/2 r^ + 24a 7 Af 3/2 r^ + 2178a 2 Af 7 r 21/2 
+ 4559a 4 AT 5 r 21/2 + 732a 6 A/ 3 r 21/2 + 81a 8 A/r 21/2 - 3336a 3 M 13/2 r 10 - 7502a 5 M 9/2 r 10 - 1328a 7 M 5/2 r 10 
+ 48a 9 VMr 10 - 786a 4 Af 6 r 19/2 + 1421a 6 Af 4 r 19/2 - 366a 8 A/ 2 r 19/2 - 1272a 3 Af 15 / 2 r 9 + 5324a 5 M n / 2 r 9 

+ 4322a 7 Af 7/2 r 9 + 174a 9 A/ 3/2 r 9 + 3501a 4 Af 7 r 17/2 - 5631a 6 Af 5 r 17/2 - 1651a 8 A/ 3 r 17/2 + 156a 10 Afr 17/2 
+ 942a 5 A/ 13 / 2 r 8 - 10a 7 A/ 9 / 2 r 8 - 1360a 9 Af 5 / 2 r 8 - 3798a 4 Af 8 r 15/2 - 5611a 6 Af 6 r 15/2 + 4900a 8 M 4 r 15/2 
+ 288a 10 M 2 rQ 5/2 + 6648a 5 M 15/2 rl - 2186a 7 A/ 11/2 r^ - 1512a 9 Af 7/2 r^ + 264a 11 A/ 3/2 r 7 1 + 4740a 6 M 7 rf /2 
+ 5292a 8 M 5 rQ 3/2 - 1686a 10 Af 3 ro 3/2 - 2088a 5 M 17/2 r% - 11734a 7 Af 13 / 2 /- 6 , + 2090a 9 M 9/2 r% 
+ 384a n Af 5/2 ^ + 2844a 6 M 8 r" /2 - 2100a 8 Af 6 r" /2 - 2490a 10 M A rf /2 + 246a 12 Af 2 r" /2 
+ 3528a 7 Af 15/2 rg + 8106a 9 A/ 11/2 r^ - 790a n Af 7/2 rg - 348a 6 M 9 rp /2 - 5398a 8 M 1 ^/ 2 + 495a 10 A/ 5 r 9/2 



o 



+ 414a 12 Af 3 r 9/2 + 360a 7 Af 17 / 2 r 4 - 2174a 9 Af 13 / 2 r 4 - 2598a n Af 9 / 2 r 4 + 120a 13 Af 5 / 2 r 4 + 708a 8 A/ 8 r 7/2 
+ 3844a 10 M 6 r 7/2 - 126a 12 Af 4 r 7/2 - 728a 9 Af 15/2 r 3 + 570a n Af 11/2 r 3 + 324a 13 A/ 7/2 r 3 - 543a 10 Af 7 r 5/2 

- 1216a 12 Af 5 r 5/2 + 24a 14 A/ 3 r 5/2 + 554a n Af 13 / 2 r 2 - 52a 13 Af 9 / 2 r 2 + 186a 12 Af 6 r 3/2 + 144a 14 Af 4 r 3/2 

- 188a 13 Af 11/2 r - 24a 14 Af 5 ^ + 24a 15 Af 9/2 ) , 

F r[2] = r l /2 ( r o /2 - 3Af^ + 2aVM^j (3Af 3 r 29/2 - 24Af 4 r 27/2 + 9a 2 A/ 2 r 27/2 - 78aAf 7 / 2 r 13 + 33A/ 5 r 25/2 

+ 57a 2 A/ 3 r 25/2 + 9a 4 A/r 25/2 + 486aM 9 / 2 r 12 + 44a 3 Af 5 / 2 r 12 + 3a 6 r 23/2 - 885a 2 AT 4 r 23/2 - 102a 4 Af 2 r 23/2 

- 612aM 11/2 r^ + 146a 3 Af 7/2 r^ + 18a 5 A/ 3/2 r^ + 1845a 2 Af 5 r 21/2 + 688a 4 M 3 r 21/2 + 9a 6 Afr 21/2 

- 1086a 3 M 9 / 2 r 10 - 354a 5 M 5 / 2 r 10 + 24a 7 VMr 10 - 1206a 2 M 6 rf' 2 - 1805a 4 Af 4 r^ 9/2 - 252a 6 M 2 r < ; 9/2 
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+ 2580a 3 M 11/2 r[J + 2546a 5 A/ 7 / 2 ^ + 66a 7 M 3/2 r^ - 528a 4 M h r 17j2 - 158a 6 Af 3 rg 7/2 + 78a 8 Air* 772 
+ 12a 3 Af 13 / 2 r 8 - 3176a 5 Af 9 / 2 r 8 - 998a 7 Af 5 / 2 r 8 - 189a 4 Af 6 ^ 572 + 3058a 6 Af 4 ^ 5/2 + 171a 8 Af 2 ^ 572 

- 24a 5 Af 11/2 r 7 + 200a 7 M 7/2 r 7 + 132a 9 M 3/2 r 7 + 1143a 4 Af 7 rf /2 - 368a 6 M 5 rQ 3/2 - 1263a 8 M 3 rf /2 

- 1962a 5 M 13/2 r 6 + 1250a 7 Af 9/2 r 6 + 216a 9 M 5/2 r 6 - 639a 6 M 6 /-" 72 - 347a 8 M 4 r" /2 + 123a 10 M 2 r" /2 
+ 684a 5 A/ 15/2 r 5 + 1978a 7 Af 11/2 r 5 - 590a 9 A/ 7/2 r 5 - 879a 6 M 7 r 9/2 + 232a 8 A/ 5 r 9/2 + 186a 10 M 3 r 9/2 

- 684a 7 M 13/2 r 4 - 832a 9 Af 9/2 r 4 + 60a n Af 5/2 r 4 + 114a 6 M 8 r 7/2 + 1048a 8 Af 6 r 7/2 - 90a 10 Af 4 r 7/2 

- 108a 7 A/ 15/2 rg + 252a 9 Af 11/2 r^ + 144a n M 7/2 rl - 147a 8 Af 7 r 5/2 - 452a 10 AfV 5/2 + 12a 12 M 3 r 5/2 
+ 146a 9 M 13/2 r 2 - 38a n M 9/2 r 2 + 69a 10 A^r 3 / 2 + 72a 12 A/ 4 r 3/2 - 70a n Af 11/2 r - 12a 12 Af 5 ^ 

+ 12a 13 M 9 / 2 ). 



E. Electromagnetic /-mode Regular izat ion Parameters 



In the electromagnetic case, the singular part of the self-force is given by 

F a =eg ah u c Af M , (2.55) 

where e is the charge of the particle, g ab is the inverse metric, u c is the four-velocity and A^) is the electromagnetic 
singular field. Here, an ambiguity arises in the definition of u a in the angular directions away from the world-line. 



In Eq. (2.551 one is free to define u a {x) as they wish provided \m\ x ^ s u a (x) = u a . A natural covariant choice would 
be to define this through parallel transport, u a (x) = g a iu b . However, in reality it is more practical in numerical 
calculations to define u a such that its components in Boyer-Lindquist coordinates are equal to the components of u a 
in Boyer-Lindquist coordinates |41j . Doing so, the regularization parameters are given by 



F t [-i] 



r r sgn Ar 

~r (a 2 + L 2 ) + 2a 2 M + r| ' 
sgn Ar (Er (a 2 + r§) + 2aM (aE - L)) 
Fr[ ~ 1] ~ (a 2 - 2Mr + r 2 ) (r (a 2 + L 2 ) + 2a 2 M + rf) ' 

F eH] = 0, F m] = 0, (2.56) 



where 

■t[0] 



F m = — ~ 2 £75 ( F m £ + F m£) ■ (2-57) 

nr 2 [r 2 + L 2 + 2 ^+ a 2 ) (2a 2 Af + a 2 r + L 2 r Q f 



Ff [Q] = -AaLM (4a 4 M 2 + 2a 4 Mr Q + 2a 2 L 2 Mr - a 2 Mr% - a 2 r 4 - L 2 r 4 ) 

+ E( - 12a 6 M 3 - 16a 6 Af 2 r - 7a 6 Air 2 - a 6 r 3 - 28a 4 L 2 M 2 r - 22a 4 L 2 Air 2 - 4a 4 L 2 r 3 - 6a 4 Af 2 r 3 
- 5a 4 Air 4 - a 4 r 5 - 15a 2 L 4 Air 2 - 5a 2 L 4 r 3 - 5a 2 L 2 Air 4 - a 2 L 2 r 5 - 2L 6 r 3 ), 
Ffi 0] = 2aLM(2a i M 2 - a 4 Mr - a 4 r 2 - a 2 L 2 Mr - 2a 2 L 2 r 2 ~ 2a 2 Air 3 - 2a 2 r 4 - L 4 r 2 - 2L 2 r 4 ) 

+ E(4a 6 M 3 + 4a 6 M 2 r + a 6 Mr 2 + 10a 4 L 2 M 2 r a + 5a 4 L 2 Mr 2 + 2a 4 M 2 r% + a 4 Mr 4 + 4a 2 L 4 Mr 2 
+ a 2 L 2 Mr 4 -a 2 L 2 4-L 4 r 5 ), 



{ F t[o] £ + F m lc ) 

? r[0] = -h^ [±J , (2-58) 

Trr 3 ( r 2 + L 2 + 2 ^- + a 2 ) (2a 2 M + a 2 r + L 2 r f (a 2 - 2Mr + r 2 ) 
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where 



F r £ [0] = L 2 (24a 6 M 4 + 28a 6 M 3 r - 6a 6 Af 2 r 2 - 11a 6 Mrl - 2a 6 r 4 + 56a 4 L 2 M 3 r + 2Aa 4 L 2 M 2 rl - \8a 4 L 2 Mrl 

- 6a 4 L 2 r 4 + 52a 4 M 3 r% + 20a 4 M 2 r 4 - 11a 4 Afr^ - 3a 4 r% + 30a 2 L 4 M 2 r 2 - 3a 2 L 4 Mr% - 6a 2 L 4 r 4 

+ A2a 2 L 2 M 2 r 4 - 5a 2 L 2 Mr% - <oa 2 L 2 r% + 8a 2 Af 2 ^ ~ 2a 2 Mr 7 - a 2 r% + 4L 6 Afr 5 - 2L 6 r 4 + 6L 4 Mr% 

- 3L 4 r% + 2L 2 Mrl - L 2 rl) 

- 2aELM(24a 6 M 3 + 36a 6 Af 2 r + 18a 6 Mr 2 + 3a 6 rg + 56a 4 L 2 M 2 r + A8a 4 L 2 Mr 2 + \Oa 4 L 2 rl + 2Aa 4 M 2 r% 
+ 24a 4 Mr 4 + 6a 4 r% + 30a 2 L 4 Mr 2 + lla 2 L 4 r% + 22a 2 L 2 Afr 4 + 9a 2 L 2 r% + 6a 2 Mr% + 3a 2 r 7 + 4L 6 rg 

+ 3£ 4 r{j + 3L 2 r 7 ) 

+ E 2 (2a 2 M + a 2 r + r%) (l2a 6 M 3 + 16a 6 Af 2 r + 7a 6 Mr 2 + a 6 ^ + 28a 4 L 2 M 2 r + 22a 4 L 2 Mrl + Aa 4 L 2 rl 
+ Qa 4 M 2 rl + ba 4 Mr 4 + a 4 4 + 15a 2 L 4 Mr 2 + 5a 2 L 4 ^ + 5a 2 L 2 Mr 4 + a 2 L 2 r% + 2L 6 r i ), 
F% 0] = -L 2 (8a 6 Af 4 + 12a 6 Af 3 r - 2a 6 Mr\ + 20a 4 L 2 M 3 r + 8a 4 L 2 M 2 rl - Aa 4 L 2 Mrl + 32a 4 M 3 ^ + 12a 4 M 2 r 4 

- Qa 4 Mrl - aV 6 , + 8a 2 L 4 M 2 rl - 2a 2 L 4 Mrl + 2Aa 2 L 2 M 2 r 4 - Aa 2 L 2 Mrl - 2a 2 L 2 r% + 8a 2 M 2 r% 

- 2a 2 Mrl ~ &2r o + 2L 4 Mr\ - L 4 r% + 2L 2 Mr 7 - L 2 r%) 

+ 2aELM(8a 6 M 3 + 12a 6 Af 2 r + 6a 6 Air 2 + a'V 3 , + 20a 4 L 2 Af 2 r + 14a 4 L 2 Afr 2 + 2a 4 L 2 r% + 16a 4 AfV 3 , 
+ 16a 4 Mr 4 + 4a 4 rg + 8a 2 L 4 M r 2 + a 2 L 4 r 3 + Ua 2 L 2 Mr 4 + 5a 2 L 2 r^ + 6a 2 Mr® + 3a 2 r 7 a + L 4 r^ + 3L 2 r 7 ) 

- E 2 (2a 2 M + a 2 r + r 3 ) (4a 6 M 3 + Aa 6 M 2 r + a 6 Mr 2 Q + 10a 4 L 2 Af 2 r + 5a 4 L 2 Afr 2 + 2a 4 Af 2 r 3 + a 4 Mr 4 
+ Aa 2 L 4 Mr 2 + a 2 L 2 Mr\ - a 2 L 2 r% - L 4 ^), 



F e[0 ] = 



(2.59) 




(2.60) 



where 



4>[0] 



Ua 4 M 2 + 11a 4 Afr + 2a 4 r 2 , + lla 2 L 2 Afr + 4a 2 L 2 r 2 l + 4a 2 AT rg 
-4a 4 Af 2 - 2a 4 Afr - 2a 2 L 2 Mr - Aa 2 Mr 3 , - a 2 ri - L 2 rt 



+ a 2 r 4 + 2L 4 r 2 +L 2 r 4 , 



As with the Scalar case, F a r 2 i proves too large to include in paper format and so is available electronically |56j . 
although we provide F r [ 2 ] for circular orbits below to allow the reader to get an understanding of the structure of the 
parameters. 



2aM 5 / 2 r 3/2 + M 2 r 2 (r - 3M) 



(2.61) 




11/2 



3 ' 



12 



where 

F r[2] = - ( r a 2 + a ^M) (-3M 3 

r o 3 ^ — 96MVq ^ — 9a 2 M 2 rg ^ + 438aM 7 / 2 rQ 5 + 495M 5 r 9 ^ — 399a 2 M 3 r 9// 

- 9a 4 Mr 29/2 - 2088aM 9 / 2 r 14 + 388a 3 M 5 / 2 r 14 - 3a 6 r 27/2 - 588M 6 r 27/2 + 1080a 2 M 4 r 27/2 
-228a 4 M 2 r 27/2 + 2058aM 11/2 r^ 3 - 372a 3 M 7/2 r 1 l 3 + 54a 5 M 3/2 r^ 3 + 4467a 2 M 5 r 25/2 + 599a 4 M 3 r 25/2 

- 123a 6 Mr 25/2 + 828aM 13/2 r^ 2 - 7504a 3 M 9/2 r^ 2 + 644a 5 Af 5/2 r^ 2 - 24a 7 VMr^ 2 - 6a 8 r 23/2 

- 11265a 2 Af 6 r 23/2 + 1578a 4 M 4 r 23/2 + 609a 6 Af 2 r 23/2 + 9576a 3 M 11/2 r 11 - 5882a 5 M 7/2 r 11 

- 528a 7 M 3/2 r 11 + 2358a 2 M 7 r 21/2 + 14389a 4 Af 5 r 21/2 + 3516a 6 M 3 r 21/2 - 189a 8 Mr 21/2 
+ 7608a 3 M 13/2 r^° - 2650a 5 M 9/2 r^° + 3296a 7 M 5/2 ^° - 48a 9 %/Mr^° - 39138a 4 M 6 r 19/2 

- 16229a 6 M 4 r 19/2 - 330a 8 Af 2 r 19/2 - 2184a 3 M 15 / 2 r 9 + 21172a 5 M n / 2 r 9 - 3530a 7 M 7 / 2 r 9 - 750a 9 M 3 / 2 r 9 
+ 21603a 4 Af 7 rJ 7/2 + 39699a 6 Af 5 r 17/2 + 8911a 8 A/ 3 r 17/2 - 156a 10 Mr 17/2 - 29982a 5 Af 13 / 2 r 8 

- 26366a 7 Af 9 / 2 r 8 + 1576a 9 M 5 / 2 r 8 - 5706a 4 M 8 r 15/2 - 16697a 6 M 6 r 15/2 - 16348a 8 M 4 r 15/2 

- 1368a 10 M 2 r 15/2 + 13800a 5 Af 15 / 2 r 7 + 39338a 7 Af n / 2 r 7 + 9096a 9 M 7 / 2 r 7 - 264a 11 Af 3 / 2 r 7 
+ 2436a 6 M 7 r 13/2 + 1368a 8 M 5 r 13/2 + 3426a 10 M 3 r 13/2 - 3096a 5 M 17/2 r^ - 23546a 7 Af 13 / 2 ^ 

- 16754a 9 Af 9/2 r 6 - 1032a 11 M 5/2 r 6 + 4404a 6 M 8 r 11/2 + 4464a 8 M 6 r 11/2 + 654a 10 M 4 r 11/2 

- 246a 12 M 2 r 11/2 + 5880a 7 M 15/2 r 5 + 15222a 9 M 11/2 r 5 + 2686a 11 M 7/2 r 5 - 516a 6 M 9 r 9/2 - 9146a 8 M 7 r 9/2 

- 3435a 10 A/ 5 r 9/2 + 54a 12 M 3 r 9/2 + 600a 7 M 17/2 r 4 - 4018a 9 Af 13/2 r 4 - 4530a 11 M 9/2 r 4 - 120a 13 M 5/2 r 4 
+ 1164a 8 M 8 r 7/2 + 6992a 10 Af 6 r 7/2 + 702a 12 AfV 7/2 - 1288a 9 M 15/2 r^ + 1158a n M 11/2 rg + 540a 13 M 7/2 r 3 1 

- 969a 10 Af 7 r 5/2 - 2336a 12 M 5 r 5/2 - 24a 14 M 3 r 5/2 + 1030a n M 13/2 r^ - 116a 13 M 9/2 r^ + 354a 12 M 6 r 3/2 
+ 288a 14 M 4 r 3/2 - 364a 13 M n / 2 r - 48a 14 M 5 ^ + 48a 15 M 9 / 2 ) , 

F r[2] = - r o ( r o /2 - SMyfro + 2a\/M) (3M 3 r 15 + 96M 4 r 14 + 9a 2 M 2 r 14 - 438aM 7/2 r 27/2 - 255Af 5 r 13 

+ 393a 2 M 3 r 13 + 9a 4 Afr 13 + 1014aAf 9/2 r 25/2 - 388a 3 M 5/2 r 25/2 + 3a 6 r 12 - 69a 2 Af 4 r 12 + 210a 4 Af 2 r 12 
+ 252aM n / 2 r 23/2 - 406a 3 M 7 / 2 r 23/2 - 54a 5 M 3 / 2 r 23/2 - 2763a 2 M 5 r 11 - 32a 4 M 3 r 11 + 105a 6 Mr 11 
+ 3594a 3 M 9 / 2 r 21/2 - 738a 5 M 5 / 2 r 21/2 + 954a 2 M 6 r 10 + 763a 4 M 4 r 10 - 204a 6 M 2 r 10 - 2940a 3 M n / 2 r 19/2 
+ 1610a 5 M 7/2 r 19/2 + 474a 7 M 3/2 r 19/2 - 3648a 4 M 5 r 9 1 - 3014a 6 A/V 9 , + 222a 8 Mr 9 , - 564a 3 M 13/2 rJ 7/2 
+ 2656a 5 M 9 / 2 r 17/2 - 2366a 7 M 5 / 2 r 17/2 + 10731a 4 A/ 6 r 8 + 12274a 6 M 4 r 8 + 963a 8 M 2 r 8 

- 21096a 5 Af n / 2 r 15/2 - 3328a 7 A/ 7 / 2 rJ 5/2 + 996a 9 M 3 / 2 r 15/2 - 1881a 4 M 7 r 7 - 7640a 6 M 5 r 7 - 9303a 8 Af 3 r 7 
+ 1 7622a 5 Af 13/2 r 13/2 + 30698a 7 M 9/2 r 13/2 + 1056a 9 Af 5/2 r 13/2 - 21807a 6 A/ 6 ^ + 517a 8 M A r% 

+ 2283a 10 A/ 2 r 6 - 1044a 5 Af 15/2 r 11/2 - 16022a 9 M 7/2 r 11/2 + 8817a 6 M 7 r 5 + 30400a 8 Af 5 r 5 + 426a 10 M 3 r 5 

- 8268a 7 M 13 / 2 r 9/2 + 4904a 9 M 9 / 2 r 9/2 + 2940a n M 5 / 2 r 9/2 - 174a 6 M 8 r 4 - 10712a 8 M 6 r 4 - 13842a 10 M 4 r 4 
+ 1428a 7 M 15/2 r 7/2 + 12492a 9 Af 11/2 r 7/2 - 360a n M 7/2 r 7/2 - 1011a 8 Af 7 r 3 + 4108a 10 M 5 r 3 + 2172a 12 M 3 r 3 

- 2062a 9 M 13 / 2 r 5/2 - 5774a n M 9 / 2 r 5/2 + 1725a 10 A/ 6 r 2 - 480a 12 M 4 r 2 + 986a n M n / 2 r 3/2 

+ 864a 13 M 7/2 r 3/2 - 876a 12 M 5 r - 17222a 7 (Mr ) 11/2 + 144a 14 M 4 - 156a 13 ^M 9 r a + 24a 7 yVr 21 ) 

F. Gravitational i-mode Regularization Parameters 

1. Self-force regularization 

The singular part of the self force on a gravitational particle is given by 

pa = k abcd h^J d , (2.62) 

where 

k abcd = l g ad u b u c _ g ab u c u d _ I^V^d + h^gbc^ + }_ g ad g bc^ (2.63) 
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and h bc is the trace-reversed singular field. Note that, as in the electromagnetic case, an ambiguity arises here due 
to the presence of terms involving the four- velocity at x. One is free to arbitrarily choose how to define this provided 
linxj.^2 u a = u a . Following Barack and Sago [31], we choose to take the Boyer-Lindquist components of the four 
velocity at x to be exactly those at x. The regularisation parameters in the gravitational case are given by, 

pt _ r f sgnAr 

ro (a 2 + L 2 ) + 2a 2 Af + rg ' { 1 

rr _ sgnAr (Er (a 2 + r 2 ) + 2aM(aE - L)) 



(a 2 - 2Mr + r 2 ) (r (a 2 + L 2 ) + 2a 2 M + rg) ' 

F» K] =0, ^[-i]=0, (2.66) 



* [°] = — ; : r§75 ^- b7 - 

^r 2 (r 2 + L 2 + + a 2j ( 2a 2 M + a 2 rQ + L 2 rQ f 

where 

F| [0 ] = AaLM (4a 4 A/ 2 + 2a 4 Air + 2a 2 L 2 Mr - a 2 Mrl - aV 4 , - L 2 r 4 ) 

+ £(l2a 6 A/ 3 + 16a 6 Af 2 r + 7a 6 Afr 2 + a 6 ^ + 28a 4 L 2 M 2 r + 22a 4 L 2 Mrl + Aa 4 L 2 rl + 6a 4 Af 2 rg + 5a 4 A/r 4 
+ a 4 rg + 15a 2 L 4 Air 2 + 5a 2 L 4 r 3 + 5a 2 L 2 Air 4 + a 2 L 2 r% + 2L e r 3 ), 
F^ [0 ] = -2aLM(2a 4 M 2 - a 4 A/r - a 4 r 2 - a 2 L 2 Mr - 2a 2 L 2 r 2 ~ 2c?Mr% - 2a 2 r\ - L 4 r 2 - 2L 2 r\) 

+ E(- 4a 6 M 3 - 4a 6 M 2 r - a 6 Air 2 - l0a 4 L 2 M 2 r - ha 4 L 2 Mrl - 2a 4 M 2 rl - a 4 Air 4 - Aa 2 L 4 Mrl 
- a 2 L 2 Mr 4 + a 2 L 2 r% + L 4 r 5 ), 



where 



[o] = : ; ^72 . ( 2 - 68 ) 

7rr 6 (r 2 + L 2 + ^M. + a 2 J (2a 2 M + a 2 r + L 2 r ) 2 (a 2 - 2A/r + r 2 ) 



i^' [0] = +i 2 rg( - 24a 6 Ai 4 - 28a 6 APV + 6a 6 Af 2 r 2 + 11a 6 Air 3 , + 2a 6 r 4 - 56a 4 L 2 Af 3 r - 24a 4 L 2 Ai 2 r 2 

+ 18a 4 L 2 Airg + 6a 4 L 2 r 4 - 52a 4 Af 3 ^ - 20a 4 M 2 r 4 + 11a 4 Mr b Q + 3a 4 r£ - 30a 2 L 4 Ai 2 r 2 + 3a 2 L 4 Mrl 
+ 6a 2 L 4 r 4 - A2a 2 L 2 M 2 r 4 + 5a 2 L 2 Mr^ + 6a 2 L 2 r% ~ 8a 2 Ai 2 r 6 + 2a 2 Mr£ + a 2 r% - AL 6 Mr 3 + 2L 6 r 4 

- QL 4 Mrl + ZL 4 rl - 2L 2 Mr 7 + L 2 rg) 
+ 2aELMr 3 (24a 6 M 3 + 36a 6 Ai 2 r + 18a 6 Mr 2 + 3a 6 r 3 + 56a 4 L 2 M 2 r + 48a 4 L 2 Air 2 + 10a 4 L 2 r 3 

+ 24a 4 Ai 2 rg + 24a 4 Air 4 + 6a 4 r^ + 30a 2 L 4 Air 2 + 11a 2 L 4 rl + 22a 2 L 2 Air 4 + 9a 2 L 2 rg + 6a 2 Air 6 , 
+ 3a 2 r 7 + 4L 6 r 3 + 3L 4 r 5 + 3L 2 r^) 

- E 2 r 3 (2a 2 M + a 2 r + r 3 ) (l2a 6 AT 3 + 16a 6 Af 2 r + 7a 6 AT r 2 + a 6 r 3 + 28a 4 L 2 M 2 r + 22a 4 L 2 Mr 2 
+ Aa 4 L 2 rl + 6a 4 A/ 2 rg + 5a 4 A/r 4 + a 4 r% + 15a 2 L 4 Mr 2 + 5a 2 L 4 rg + 5a 2 L 2 A/r 4 + a 2 L 2 rl + 2L 6 

F£ [0] = -L 2 rl( - 8a 6 M 4 - 12a 6 M 3 r + 2a 6 Mr 3 , ~ 20a 4 L 2 M 3 r - 8a 4 L 2 M 2 rl + Aa 4 L 2 Mrl - 32a 4 AfV 3 , 

- 12a 4 A/ 2 r 4 + 6a 4 A/r^ + a 4 rg - 8a 2 L 4 A/ 2 r 2 + 2a 2 i 4 Afrg - 24a 2 L 2 A/ 2 r 4 + 4a 2 i 2 Mr^ + 2a 2 L 2 r^ 

- 8a 2 Af 2 r 6 + 2a 2 AfrJ + a 2 r« - 2L 4 A7r^ + L 4 r% - 2L 2 Mr 7 + L 2 r*) 

- 2aELMrl(8a & M 3 + 12a 6 M 2 r a + 6a 6 Air 2 + a 6 r^ + 20a 4 £ 2 A/ 2 r + 14a 4 L 2 Afr 2 + 2a 4 L 2 r 3 l + 16a 4 M 2 rg 
+ 16a 4 Mr 4 + 4a 4 r^ + 8a 2 L 4 Mrl + a 2 L 4 rl + Ua 2 L 2 Mr 4 + ha 2 L 2 r% + 6a 2 Mr% + 3a 2 r 7 + L 4 r% + 3L 2 r 7 ) 

- E 2 r 3 (2a 2 M + a 2 r + rg) ( - 4a 6 A/ 3 - 4a 6 Af 2 r - a 6 A/r 2 - 10a 4 L 2 A/ 2 r - 5a 4 i 2 A/r 2 - 2a 4 A/ 2 r^ 

- a 4 Air 4 - 4a 2 i 4 Mr 2 - a 2 L 2 Afr 4 + a 2 L 2 rl + L 4 ^), 



r6 r 3 ), 



F [o]=0 



(2.69) 
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F V] - — ; ^ " " • ^ 2 - 7 °) 

7rr 6 (r 2 + L 2 + + a 2 J (2a 2 AT + a 2 r + L 2 r ) 2 (a 2 - 2Mr + r 2 ) 
where 

F e[0] = -14a 4 Af 2 - lla 4 Afr - 2a 4 r 2 - lla 2 L 2 Afr - 4a 2 L 2 r 2 - 4a 2 A/r 3 - a 2 r 4 - 2L 4 r 2 - L 2 r 4 , 

F$ [0] = 4a 4 M 2 + 2a 4 Mr + 2a 2 L 2 Mr + 4a 2 Mr% + a 2 r 4 + L 2 r 4 . 

As with the scalar and electromagnetic cases, is too large for paper format and so is available electronically |56| . 
Instead, we give here FT 2 -\ for circular orbits, 



b l 2 l - ^172 



2aM 5 / 2 rn /2 + M 2 r 2 (r -3Af) r ,— OJlr J- 2 

2aVM + y^(r - 3M) 



^ (a\/M + r 3/2 ) [a 2 + r (r -2Af)] 

[a 2 + r (r - 2Af)]^ (fj [2] £ + F£ [2] /C) 

(a 4 A/ + 2a 3 v / Mro + a 2 r (-2Af 2 + Mr + r 2 ) - 4aM 3 / 2 r 5/2 + ATr 4 ) 



2 ' 



where 
F, 



£ [2] = - (a\fTl + r 3 / 2 ) (a 2 - 2Mr + r 2 ) (32a 11 Af 5/2 - 48a 10 Af 3 ^ + 144a 10 Af 2 r 3/2 + 478a 9 A/ 3/2 r 3 

- 104a 9 Af 5/2 r 2 - 192a 9 Af 7/2 r + 288a 8 AfV 3 / 2 - 992a 8 A'/ 3 r 5/2 - 473a 8 Af 2 r 7/2 + 779a 8 Mr 9/2 

- 292a 7 Af 3/2 r 5 , - 2496a 7 A/ 5/2 r 4 + 800a 7 Af 7/2 r^ + 384a 7 A/ 9/2 r 2 + 540a 7 VMr% - 576a 6 Af 5 r 5/2 
+ 2240a 6 A/ 4 r 7/2 + 2489a 6 AfV 9/2 - 4623a 6 A/V" 72 + 227a 6 Afr" /2 + 127a 6 r 45/2 - 1846a 5 M 3/2 rl 
+ 3574a 5 Af 5 / 2 ^ + 4018a 5 M ll2 r\ - l952a 5 M 9/2 r 4 - 256a 5 M 11/2 r^ + 386a 5 VMr% + 384a 4 M 6 r 7 /2 

- 1664a 4 Af 5 r 9/2 - 3950a 4 Af 4 r 41/2 + 6038a 4 M 3 r 1 ^ 2 - 2275a 4 Af 2 ^ 5/2 - 462a 4 Mr 47/2 + 105a 4 r 49/2 

- 14a 3 AT 3/2 r[J + 4444a 3 M 5/2 r% - 65Ua 3 M 7/2 r 7 - 1876a 3 M 9/2 r% + 1536a 3 M lx ' 2 r% - 144a 3 %/Mr 1 , 
+ 1728a 2 M 5 rf /2 + 674a 2 A/ 4 r 45/2 + 588a 2 Af 3 ^ 772 - 1312a 2 A/ 2 r 49/2 + 348a 2 A/r 21/2 - 22a 2 r 23/2 

- 148aAf 3/2 r^ + 884aM 5/2 r 1 1 ° - l370aM 7/2 r% + 384aA/ 9/2 r 8 ) - 192Af 4 ro 9/2 + 287M 3 r 21/2 

- 131A/ 2 r 23/2 + 18Mr 25/2 ), 



F£ [2] = r 3 , [2aVM - 3My/r^ + r 3/2 j (l44a 12 A/ 2 - 96a 11 y/M 5 r + 576a n (Mr ) 3/2 - 921a 10 A/ 3 r + 81a 10 M 2 r 2 

+ 864a 10 A/r 3 + 862a 9 A/ 3/2 r 7/2 + 582a 9 Af 7/2 r 3/2 + 576a 9 ^A/r 9 - 4010a 9 (Mr ) 5/2 + 2052a 8 AfV 2 
+ 313a 8 AfV 3 , - 5342a 8 M 2 r 4 + 1351a 8 Mr^ + 144a 8 r 6 l - 2308a 7 A/ 3/2 rQ 1/2 - 3816a 7 A/ 5/2 r 9/2 

- 1176a 7 A/ 9/2 r 5/2 + 920a 7 Mr™ + 9524a 7 (Air ) 7/2 - 1764a 6 AT 5 ?- 3 , - 2260a 6 Af 4 r 4 + 12726a 6 AfV 5 , 

- 6148a 6 Af 2 r 6 + 60a 6 A/r 7 + 254a 6 r 8 - 2580a 5 M 3/2 r 15/2 + 4204a 5 Af 5/2 r 13/2 + 3504a 5 Af 7/2 r 11/2 



+ 792a 5 M 11/2 r 7/2 + 204a E y A/r 17 - 8384a 5 (Afr ) 9/2 + 384a 4 Af 6 r 4 + 2620a 4 A/V 5 , - 11878a 4 AfV 6 
+ 8425a 4 AfV 7 , - 1656a 4 AfV 8 , - 111a 4 Afr 9 , + 88a 4 rJ° - 242a 3 Af 3/ V 19/2 + 2928a 3 Af 5/ V 17/2 



4866a 3 M 7,2 rf /2 + 1372a 3 M 9,2 rf /2 - lUa 3 \j Mr 21 + 1536a 3 (A/r ) 11/2 + 1728a 2 AfV 7 , - 1614a 2 M 4 rg 

1373a 2 M 3 rl - 1075a 2 Af 2 r^° + 334a 2 Afr 1 , 1 - 22a 2 r 1 , 2 - 148aA/ 3/2 r 23/2 + 718aAf 5/2 r 21/2 
1040aAf 7/2 rg 9/2 + 384aAf 9/2 r 47/2 - 192Af 4 r^° + 280M 3 r^ - 128M 2 r^ 2 + 18Afr^ 3 ). 



2. huu regularization 

The quantity 



H {R) = i*2W (2.71) 
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was first proposed by Detweiler [57] as a tool for constructing gauge invariant measurements from self-force calculations. 
It has since proven invaluable in extracting gauge invariant results from gauge dependent self- force calculations [39|I58|. 

Much the same as with self-force calculations, the calculation of requires the subtraction of the appropriate 

singular piece, u a u from the full retarded field. In this section, we give this subtraction in the form of 

mode-sum regularization parameters. In doing so, we keep with our convention that the term proportional to I + \ 
is denoted by (= in this case), the constant term is denoted by -ff[o], and so on. 

Note that, as in the self- force case, an ambiguity arises here due to the presence of terms involving the four- velocity 
at x. One is free to arbitrarily choose how to define this provided lim^^j u a — u a . As before, we choose this in 
such a way that the Boyer-Lindquist components of the four velocity at x are exactly those at x. The regularisation 
parameters are then given by, 



where 



H [0] = — — (2.72) 



H m = 0, (2.73) 



hLs + h&k 

H m z v " .J ' ' ( 2 - 74 ) 

Z-ktI (r 2 + L 2 + + a 2 j ( 2a 2 M + a 2 rQ + L 2 rQ f 



H[ 2] = (l2Mr^a 12 + 92M 2 r 4 a 12 + 264Af 3 r^a 12 + 336Af 4 r 2 a 12 + 160Af 5 r a 12 - 24rga 10 - 2A0Mr 7 a 
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- 1104L 2 M 6 a 10 - 1104M 2 r> 10 - 2880M 3 rfja 10 + 48L 2 Mr^a w - 4272Af 4 r 4 a 10 + 230L 2 M 2 r 4 a 10 

- 3264A/ 5 r 3 ) a 10 + 96i 2 M 3 r^a 10 - 960Af 6 r^a 10 - 1116L 2 M 4 r 2 a 10 - 2096L 2 M 5 r a 10 - 48r 1 ,V 

- 420Af r^a 8 - 120L 2 r^a 8 - 1556M 2 r^a 8 - 2872A/ 3 r£a 8 - 882L 2 Mr 7 a s - 2448A/ 4 r> 8 - 2781L 2 Af 2 r£a 8 

- 672M 5 t> 8 - 4770L 2 Af 3 r*a 8 + 72L 4 Mr^a 8 - 4272L 2 A/ 4 r 4 a 8 + 90L 4 AT 2 r 4 a 8 - 1440L 2 Af 5 r 3 a 8 

- 1044 J L 4 Af 3 r 3 1 a 8 - 2928L 4 Af 4 r^a 8 - 2112L 4 AT 5 r a 8 - 24r 1 , 2 a 6 - 168A/r^a 6 - 195L 2 rJ a 6 - 456A/ 2 r 1 ,°a 6 

- 480A/ 3 r^a 6 - 1086L 2 A/r^a 6 - 240£ 4 r 8 ,a 6 - 96A/ 4 r^a 6 - 2119L 2 Af 2 rga 6 - l528L 2 M 3 r 7 a 6 

- 1098L 4 Af r 7 a 6 - 84L 2 M 4 r 6 1 a 6 - 1578 J L 4 A/ 2 r^a 6 - 696L 4 Af 3 r^a 6 + 48 J L 6 Air^a 6 + 84L 4 Af 4 r 4 a 6 

- 190L 6 Af 2 r 4 a 6 - 1320L 6 M 3 r%a 6 - 1476L 6 A/ 4 r 2 a 6 - 75L 2 r^ 2 a 4 - 246L 2 Air 1 , 1 a 4 - 297 J L 4 r 1 ,°a 4 

- 84L 2 A/ 2 r 40 a 4 + 216L 2 Af 3 r^a 4 - 690L 4 Air^a 4 - 240L 6 r 8 a 4 + 529L 4 A/ 2 r 8 a 4 + 1374L 4 Af 3 r^a 4 

- 402L 6 Air V + 771L 6 M 2 r 6 ,a 4 + 1194L 6 M 3 r^a 4 + U^Mr^ - 190L 8 Af 2 r^a 4 - 444L 8 Af 3 rfja 4 

- 78L 4 r 1 , 2 a 2 + ZQ^Mr^a 2 - 201L 6 r 1 1 °a 2 + 384L 4 Af 2 r^°a 2 + 198L 6 Af r^a 2 - 120L 8 r§a 2 + 1092L 6 Af 2 rga 2 
+ 162L s Mr 7 a 2 + 672L 8 M 2 r 6 ,a 2 - 48L 10 AT 2 r 4 a 2 - 27L 6 r 1 , 2 + 114L 6 A/r^ - 51L 8 r^° + 222L s Mrl 

- 2AL w rl + 108i 10 Afr^) 

AaELM ( - 71a 4 r 1 , 1 - 75L 4 r 1 , 1 - 146a 2 J L 2 r^ 1 - 262a 4 A/r 1 , - 270a 2 L 2 Air 1 , - 147a 6 r[J - 159L 6 r^ 

- 465a 2 L 4 r^ - 453a 4 L 2 rg - 240a 4 Af 2 ^ - 811a 6 Air 8 , - 873a 2 L 4 Air 8 , - 1684a 4 L 2 Air 8 ) - 76a 8 r£ - 78L s r 7 

- 310a 2 L 6 rl - 462a 4 L 4 r^ - 306a 6 L 2 r\ - 1490a 6 M 2 r£ - 1542a 4 L 2 A/ 2 r 7 - 912a 6 A/V 6 , - 522a 8 Mr% 

- 549a 2 L 6 Alr 6 , - 1620a 4 i 4 Afr^ - 1 593a 6 L 2 Air 6 - 1323a 8 Af 2 r^ - 1368a 4 L 4 Af 2 r% - 2691a 6 L 2 M 2 r 3 

- 1460a 8 M 3 r 4 - 1458a 6 L 2 Af 3 r 4 + 37a 10 Air 4 + 24a 2 L 8 Air 4 + 109a 4 L 6 Afr 4 + 183a 6 L 4 Air 4 

+ 135a 8 L 2 Air 4 - 588a 8 AfV 3 , + 291a 10 Af 2 rg + 222a 4 L 6 M 2 r 3 , + 735a 6 L 4 A/ 2 r 3 , + 804a 8 L 2 M 2 r^ 
+ 858a 10 Af'V 2 + 738a 6 L 4 Af 3 r 2 + 1596a 8 L 2 Af 3 r 2 + 1124a 10 Af 4 r + 1056a 8 L 2 M 4 r Q + 552a 10 M 5 ) 
3^ 2 (368A/ 6 a 12 + 4Air 3 a 12 + 55Af 2 r 4 a 12 + 280Af 3 r^a 12 + 680M 4 r 2 a 12 + 800A/ 5 r o a 12 - 9rga 10 

- 124Afr^a 10 - 604A/ 2 r 6 a 10 - 1368M 3 rj>a 10 + 12 J L 2 Afrga 10 - 1472M 4 r 4 a 10 + 160L 2 Af 2 r 4 a 10 

- 608A/ 5 r^a 10 + 672i 2 Af 3 rga 10 + 1152L 2 M 4 r 2 a 10 + 704L 2 Af 5 r a 10 - 18r^°a 8 - 224Mrga 8 - 35L 2 rga 8 

- 901A/ 2 rga 8 - 1492Af 3 r^a 8 - 438L 2 Mr^a 8 - 884Af 4 r^a 8 - 1685L 2 Af 2 rga 8 - 2604 J L 2 AT 3 rga 8 

+ 12L 4 Afr^a 8 - 1412L 2 M 4 r 4 a 8 + 171i 4 Af 2 r 4 a 8 + 540L 4 M 3 rga 8 + 492L 4 AT 4 r 2 a 8 - 9r^ 2 a 6 - 96M r^a 6 
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- 55L 2 rl°a 6 - 278M 2 r 1 ,°a 6 - 2UM 3 r^a 6 - 596L 2 Mr^a 6 - 51L 4 rga 6 - 1679L 2 M 2 r 8 ) a 6 - IA1AL 2 M s r 7 a 6 

- 572L 4 Mr 7 a 6 - 1567L 4 Af 2 r> 6 - 1254L 4 A/ 3 r^a 6 + 4L 6 Mr%a 6 + 82L 6 M 2 r 4 a 6 + 148Z e M 3 rga 6 

- 20L 2 r 1 1 V - UG^Mr^a 4 - 56L 4 rl°a 4 - 266L 2 M^a 4 - 512L 4 Mr 9 a 4 - 33L 6 r«a 4 - 774L 4 M 2 r*a 4 

- 326L e Mr 7 a 4 - 486L 6 M 2 r%a 4 + 16L 8 M 2 r 4 a 4 - 13£ 4 r£ 2 a 2 - 80L 4 MrJV - 19L 6 rJV - UQL 6 Mr%a 2 

- 8L 8 r%a 2 - 68L 8 Mr 7 a 2 - 2L 6 r^ 2 ), 

Hfa = rg(24r^a 10 + 180Af r 4 a 10 + 512M 2 rj|a 10 + 656AT 3 r 2 a 10 + 320Af 4 r a 10 + 93r 7 a 8 + 624A/ t> 8 - 912i 2 M 5 o 8 
+ 120i 2 r^a 8 + 1448A/ 2 rga 8 + 976A/ 3 r 4 a 8 + 636L 2 Alr 4 a 8 - 912Af 4 r^a 8 + 891L 2 Af 2 r^a 8 - 1152M 5 r 2 a 8 

- 434£ 2 A/ 3 r 2 a 8 - 1720L 2 Af 4 r a 8 + 69t> 6 + 354A/r^a 6 + 375L 2 r 7 a 6 + 492Af 2 r^a 6 - 168A/ 3 r^a 6 
+ 1620L 2 Mr^a 6 + 240L 4 ^a 6 - 576M 4 rga 6 + 1207L 2 Af 2 r^a 6 - 2720L 2 M 3 r 4 a 6 + 732L 4 M r 4 a 6 

- 3372L 2 M 4 r 3 1 a 6 - 450L 4 A/ 2 r 3 1 a 6 - 2896L 4 Af 3 r 2 a 6 - 2052L 4 A/ 4 r a 6 + 21CLL 2 ^a 4 + 534L 2 Af rga 4 
+ 567L 4 r£a 4 - 384L 2 Af 2 r^a 4 - 1224L 2 A/ 3 r> 4 + 1074L 4 Mr^a 4 + 240L 6 r^a 4 - 2017L 4 Af 2 rga 4 

- 3726L 4 Af 3 r 4 a 4 + 180L 6 A/r 4 o 4 - 1525L 6 M 2 rga 4 - 1806£ 6 M 3 r 2 a 4 + 213L 4 r[Ja 2 - 18L 4 Afr^a 2 
+ 381L 6 r^a 2 - 888L 4 M 2 r 7 a 2 - 216L 6 Mr%a 2 + 120L 8 r^a 2 - 1776L 6 M 2 r^a 2 - 192L 8 A/r 4 a 2 

- 696L 8 M 2 r^a 2 + 72L 6 r% - 198L 6 Air 8 + 96L 8 r 7 - 294L 8 Af + 24L 10 rfj - 96L 10 Af r 4 ) 

+ \ar\ELM (456AT 4 a 8 + 45r 4 a 8 + 322A/r 3 ) a 8 + 862Af 2 r 2 a 8 + 1024M 3 r a 8 + 116r^a 6 + 636AT r%a 6 

+ 183L 2 r 4 a 6 + 1162Af 2 r 4 a 6 + 708A/ 3 r 3 ,a 6 + 992L 2 Mr 3 1 a 6 + 1765L 2 A/ 2 r 2 a 6 + 1026L 2 Af 3 r a 6 + 71r 8 ,a 4 
+ 262A/r£a 4 + 358L 2 rga 4 + 240M 2 r^a 4 + 1326L 2 Mrjja 4 + 279L 4 r 4 a 4 + 1206L 2 Af 2 r 4 a 4 + 10 18L 4 Mr 3 , a 4 
+ 903L 4 A/ 2 r 2 a 4 + 146L 2 r 8 a 2 + 270L 2 Mr 7 a 2 + 368£ 4 r^a 2 + 690L 4 A/r[ja 2 + 189L 6 r 4 a 2 + 348L 6 Af r 3 a 2 
+ 7bL 4 rl + l26L 6 r% + 48L 8 r 4 ) 

- 3r 3 £ 2 (304AfV + 8r^a w + 90Mr 4 a 10 + 386A/ 2 r 3 a 10 + 796A/ 3 r 2 a 10 + 792AT 4 r a 10 + 32r£a 8 

+ 308Afr^a 8 + 32L 2 r^a 8 + 1068M 2 r^a 8 + 1600A/ 3 r 4 a 8 + 304L 2 Afr 4 a 8 + 880A/ 4 r 3 a 8 + 1003L 2 Af 2 r 3 a 8 
+ 1388L 2 Ar 3 r 2 a 8 + 684L 2 M 4 r a 8 + 24^a 6 + 186Afr 8 a 6 + 113.L 2 r£a 6 + 458Af 2 r^a 6 + 364A/ 3 r°a 6 
+ 854L 2 Af r^a 6 + 48L 4 r^a 6 + 2017L 2 Af 2 r^a 6 + 1522 J L 2 Af 3 r 4 a 6 + 370L 4 Mr 4 a 6 + 849L 4 Af 2 r^a 6 
+ 602£ 4 A/ 3 r 2 a 6 + 65L 2 ^a 4 + 356L 2 Af rga 4 + 146L 4 r£a 4 + 446L 2 A/ 2 r^a 4 + 77%L 4 Mr%a 4 + 32L 6 r%a 4 
+ 945L 4 M 2 r%a 4 + l88L 6 Mr 4 a 4 + 232L 6 M 2 r 3 ] a 4 + 58L 4 r> 2 + 170L 4 Af rga 2 + 81L 6 ^a 2 + 230L 6 Mr%a 2 
+ 8L 8 rla 2 + 32L 8 Mr 4 a 2 + l7L 6 r% + 16L 8 r 7 ). 

G. Example 

To illustrate the effectiveness of the higher-order regularization parameters, we consider as an example the case of a 
scalar charge on an eccentric geodesic orbit with E — 0.955492 and L = 3.59656 in a Kerr spacetime with a = 0.5M. 
The self-force in this case was computed in Ref. [52] using a frequency domain calculation of the retarded field in 
combination with the first two regularization parameters. Figure [l] shows the effect of using higher order regularization 
parameters on this calculation. As expected, the numerical ^-modes computed in Ref. [45] asymptotically fall off as 
l~ 2 . Our FL2] regularization parameter analytically gives the coefficient of this leading order in \jl behaviour. After 
subtracting this leading order behaviour from the numerical modes, we find a remainder which falls off as l~ 4 , as 
expected, with the coefficient given analytically by our regularization parameter. Upon subtraction of F l r ^ the 

remainder falls off as l~ 6 , as anticipated. 

III. EFFECTIVE SOURCE AND m-MODE REGULARIZATION 
A. Effective Source Approach to the Self-force 

The effective source approach - independently proposed by Barack and Golbourn [15] and by Vega and Detweiler 
[16] - relies on knowledge of the singular field to derive an equation for a regularized field that gives the self-force 
without any need for post-processed regularization. If the singular field is known exactly, then the regularized field 
is totally regular and is a solution of the homogeneous wave equation. In reality, exact expressions for the singular 
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1-modes of the Regularized Scalar Self-force 




2 3 4 5 7 9 11 14 18 



FIG. 1. Regularization of the i-modes of the radial component of the scalar self- force in Kerr spacetime for the case of a 
particle in an eccentric, equatorial geodesic orbit with E = 0.955492 and L — 3.59656. The background Kerr black hole has 
spin a — 0.5M. In decreasing slope the lines represent the unregularized self-force (yellow), and the self-force regularized by 
subtracting from it in turn the cumulative sum of ^Vf-i] (blue), iOo] (green), FL2] (red) and i^[4] (black). 

field can only be obtained for very simple spacetimes. More generally, the best one can do is an approximation such 
as the one derived in this paper. Splitting the retarded field into approximate singular and regularized parts, 



^fret) = £(S)+£(R): 



(3.1) 



substituting into the wave equation, Eq. (2.1) of Paper I, and rearranging, we obtain an equation for the regularized 
field, 



T>A ~B qA 



with an effective source, 



S. 



off 



-V 



A Bvf S) ~^Q I U A 6 4 (x, Z (T'))dT'. 



(3.2) 



(3.3) 



For sufficiently good approximations to the singular field, and S A arc finite everywhere, in particular on the 

world-line. As a result, one never encounters problematic singularities or <5-functions, making the approach particularly 
suitable for use in time domain numerical simulations. A detailed review of this approach can be found in |591 160j . 

One disadvantage of the effective source approach stems from the fact that the source must be evaluated in an 
extended region around the world-line. Since the source is derived from a complicated expansion approximating 
the singular field, its evaluation can dominate the run time of a numerical code. This problem is exasperated as 
increasingly good approximations to the singular field — using increasingly high order series expansions — are used, 
placing a practical upper limit on the order of singular field approximation which may be used in effective source 
calculations. Existing calculations [61 63 settled on what appears to be a "sweet spot", using an approximation 
accurate to 0(e 2 ). 

Despite it being possible to compute higher order effective sources from our singular field approximation this 
practical consideration may appear to rule out the usefulness of high-order expansions of the singular field in effective 
source calculations, particularly in the case of the Kerr spacetime where even an order 0{e 2 ) approximation to the 
singular field is quite unwieldy. However, it turns out that high-order expansions can, in fact, be put to good use in 
effective source calculations. In this section, we show how this may be achieved in the case of the m-mode approach 
to effective source calculations. In this approach, one first performs a decomposition into m-modes 



1 

2^ 



-im4> 



d<t>. 



(3.4) 



and independently evolves the m-decomposed form of the wave equation for each m-mode. These equations have an 
m dependent effective source which is derived from the particular choice of approximation to the singular field. The 
full field is then given as a sum of these individual modes: 



$ = 



m=— o 



(3.5) 
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For an approximation accurate to 0(e n ), the numerical solutions for the field fall off as m,~( n+2 ) for m even and as 
m -(n+3) f or m ^ Obviously, only finitely many m- modes (typically ~ 10 — 20) can ever be computed numerically; 
with the error from truncating the sum at a finite m putting an upper limit on the accuracy of the self-force which 
can be computed. This may be mitigated somewhat by fitting for a large-m tail, but that fit itself requires more 
modes and is only ever approximate. Here, we propose a much better solution, using the higher-order terms in the 
singular field (those which have not been used in computing the effective source) to analytically deriving expressions 
for the tail. In many ways, this is analogous to the /-mode regularization scheme where there is a large-/ tail and one 
can compute /-mode regularization parameters. 



B. Derivation of m-mode Regularization Parameters 

To derive analytic expressions for the large-m tail, we first note that an approximation to the singular field accurate 
to 0(e n ) can be written in the form 



_ 3(71+1) 



3(71 + 1) 



^H-sUf E A l ,7sin l (A0/2)+ A n , iS in i - 1 (A0/2) S in(A^)l +0(e n+1 ) 

Po i=o i=o 



i=0 
.3(71+1) 



3(71+1) 



2n 

Po 



J2 ^n,7sin l (A0/2) + 2 J2 ^sin i; (A</>/2)cos(A0/2)] +0(e n+1 ) 

i=0 i=0 

odd 



(3.6) 



where the coefficients A„ ti are functions of position ro and 9 0l of the constants of motion E, L, C, and also of Ar 
and A8. This form has the benefit of ensuring that the approximation is regular everywhere except on the world-line, 
while still being amenable to analytic integration in the <f> direction. This makes it particularly appropriate for use in 
m-mode effective source calculations [64] . 

Using the leading orders (say, to 0(e p )) in this expansion to compute an effective source, one is left with a singular 
field remainder which is finite, but of limited differentiability on the world-line. Since it is finite, we can safely set 
Ar = A9 = in Eq. (3.6), leading to a singular field remainder which has the form 



<f> s (x) 



n n 

20(A0)-ll[ 5 n , i sin i (A^/2) + 2 ^ B ^ shi ,; (A0/2) cos(A0/2) +C(e" +1 ), (3.7) 



where 0(A0) is the Heaviside step function. Substituting this into Eq. (3.4) and noting that for even j 

2im 



20(A<£) - 1 sin J '(A<£/2)cos(A0/2)e _im *# = 



3 + 1 J-* 



20(A0)-1 sin J+1 (A0/2) e - mi0 d0 (3. 



we are left with trivial integrals of the form 



20(A0) - 1 



J+ 1 



(A(/>/2)e- lm< ^ : 



20(A0)-1 sin 3+1 (A^/2)cos(m- 



(3.9) 



As a result, we see that the real regularization parameters are given by the odd terms in the expansion of the singular 
field and the imaginary parameters are given by the even terms. Furthermore, we see that the falloff with m is always 
an even power of 1/ro in the real part and an odd power of 1/m in the imaginary part. 

While this analysis was done for the field, it should be noted that it equally well applies to the self-force. The only 
modification necessary is to compute the self-force from the singular field before setting Ar = A9 = 0; the remainder 
of the calculation proceeds in exactly the same way. 

Finally, we note the m-mode regularization parameters derived in this way are dependent on the singular field being 
written in the form given in Eq. (3.6). Effective source calculation may use some other form for the approximation 



to the singular field (while still being accurate to the same order), in which case there is no guarantee that the 
regularization parameters given here are appropriate. 
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C. m-mode Regularization Parameters 

Below, we give the results of applying this calculation to each of scalar, electromagnetic and gravitational cases in 
turn. In doing so, we omit the explicit dependence on m which in each case is 

-4-F a [2] „ m 24F a[4] 



a[2] ~ 7r(2m - l)(2m + 1) ' a[4] ~ 7r(2m - 3) (2m - l)(2m + l)(2m + 3) ' 

pm = -480F n[6] 

a[6] tt(2to - 5) (2m - 3) (2m - l)(2m + l)(2m + 3) (2m + 5) ' 



20160F a 

7r(2m - 7) (2m - 5) (2m - 3) (2m- l)(2m + l)(2m + 3) (2m + 5) (2m + 7) ' 



Kk = _,o„ ,n/o„ ^ on,o„ ^muo-.o^,^,^ ( 3 - 10 ) 



As the expressions for generic orbits of the Kerr spacetime are too large to be of use in printed form we give here 
only the representative expressions for two cases: the r component of the scalar self-force for a circular geodesic orbit 
and the quantity H — ^h a i,u a u b in the gravitational, circular orbits case. We direct the reader online |65j for more 
generic expressions in electronic form. 

For circular orbits, the scalar self-force m-mode regularization parameters are: 



F = M 

r[2] 24r^[aM + ro^/Mra^a^/MFo + r {r - 3M)] 3 / 2 [« 2 + r (r - 2M)] 3 / 



24a 7 M 2 



2 

24a 6 M v / Mr (M - 2r ) - 4a 5 Mr (23M 2 + Mr - 6r%) + 2a 4 Mr v / Mr (45M 2 - 112Mr + 31r 2 ) 

'o; 
(3.H) 



2a 3 Mr 2 (45M 3 + 45M 2 r - 73Mr 2 + 19r 3 ) - 3a 2 r 2 y^M^(29M 4 - 88M 3 r + 38M 2 r 2 - 4Mr 3 + 
6aMr 4 (29M 3 - 43M 2 r + 21Mr 2 - 3r 3 ) - 3r|j y/Mr (29M 3 - 25M 2 r + 3Afr 2 + r 3 ) 



F, 



M 2 



r[4] 



23040a 14 M 2 v /]Ur 



1440rg[aA/ + r VMro] [W^ro + r (r - 3M)] 7 / 2 [a 2 + rQ ( rQ _ 2M)] 3 / 2 
+ 11520a 13 Af 2 r (Af - 8r ) + 384a 12 Afr v /M7 ; (461Af 2 - 81Afr - 361 ,, 
- 3780A/ 2 r + 1233Afr 2 + 480r 3 ) - 64a 10 r 2 ^/m^(8M9M 4 - 3593Af 3 r - 16212Af 2 r 2 
+ 6336Mrg + 360r 4 ) + 32a 9 Mr 3 ) (2835M 4 - 69401A/ 3 r + 46565Af 2 r 2 + 13779 AfrfJ - 8748r, 
+ 192a s rly/Mro(U70M 5 - 3621Af 4 r - 15645M 3 
+ 16a 7 r 4 (-1479A/ 6 + 210966M 5 r - 224760M 4 r 2 
16a 6 r 4 VA/r (43101A/ 6 - 61443M 5 r 

>■■■■■>> ■ 337457M 5 r 2 + 71894M 4 r^ - 262111Af 3 r 4 + 111498Af 2 ^ 



u Mr 2 (307M 3 



12662Af 2 r 3 l 
- 49213M 3 r 3 l 



+ 4158rg) + 12a 5 r^(-2367Af 7 - 221220M 6 r 4 

1 AA ecn )l/f»,6 i c;qq„7\ , o„4„5 . / ji..r„_ rtRi ok Jl/f7 



1529Afr 4 - 342r^) 
93619M 2 r 4 - 19953Mr^ + 180r£) 



+ 



343776A/ 3 r 3 - 100489A/ 2 r 4 - 7763Mr 



■ 1949709M 4 rjj 



- 14459Afrg + 588r 7 ) + 3a 4 r^ v/M^(76125A/ 7 - 176307Af 6 r - 1157559M 5 r 2 + 19497U9A, , 

- 855873Af 3 r 4 - 26505M 2 r^ + 76235Mrg - 806br 7 ) + 12a 3 r^(76125Af 7 - 152637M 6 r - 93174M 5 r, 
+ 281414Af 4 rg - 166063Af 3 r 4 + 36555M 2 rfj - 3480Afr^ + A60r 7 ) + I8a 2 4 v / Aini(76125Af 6 

1 A C1 oo 71 ,f 5„ i 711771 7i,f4 2 i i cGfic )\/r3„3 i nnnc Jl ,f 2 4 , oinnuJ i i<">c„6\ i net — 10/oc 



"5„2 



16696M 3 rg 



- 3190Mr£ 



- 145182A/ 5 r + 70771Af 4 r^ 

- 47369Af 5 r + 31856A/ 4 r 2 - 8692Af 3 rg + 705M 2 r 4 - 75Afrfj + 40r^) 



19905M 2 r 4 



-225r^) + 36ar(i (25375Af 6 



9r^ 1 v /A77 ; o(25375Af 5 - 47015M 4 r + 29014A/ 3 r 2 - 4814Af 2 rg - 1365Afr 4 + 405r^) . (3.12) 



and the gravitational, m-mode parameters for H are 



H[2] = 



\ 



M (a 2 +r (r -2M)) 



2aVM + ^/7 ; o(r - 3M) 

44a 4 M + 88a 3 y/M^ - 3a 2 r (M - r )(29Af + 15r ) + 6ayMro /2 (14r - 29M) - 87Mr^ + 45r§ 



12r, 



15/4 



(a 2 + r (r - 2M)) 



(3.13) 
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m modes of the radial self-force 




m 



FIG. 2. Regularization of the radial component of the self-force for the case of a scalar particle on a circular geodesic of 
radius ro = 10M in Kerr spacetime with a = 0.6M. The numerical self-force modes (blue dots) asymptotically match the Fmi 
regularization parameter (blue line) for large m. After regularization, the remainder (purple squares) fall off as m -6 . 



D. Example 

As an example application of these m-mode regularization parameters, we consider the case of a scalar charge on a 
circular geodesic orbit of radius 10M in the Kerr spacetime with a = 0.6M. The self-force in this case was computed 
in Ref. |64j using the m-mode effective source approach with an effective source derived from an approximation to 



the singular field of the form (3.6) accurate to 0(e 2 ). As expected, this gave numerical results for the m-modes of 
the self-force which asymptotically fall off as mT A . In this case, the i 7 ^, parameter is not needed as it has already 
been subtracted through the effective source calculation. However, the parameter has not been subtracted and 
asymptotically gives the leading order behaviour (in 1/m) of the modes. Subtracting this from the numerical results 
therefore leave a remainder which falls off as m~ 6 . Furthermore, a numerical fit of this remainder can be done to 
numerically determine the next two parameters, in this case giving F r rgi = 0.108797 and F r ig\ = 11.3398. 

In Fig. |2j we plot the results of subtracting the i 7 ^ and F™ 6 j regularization parameters in turn from the raw 

numerical data. For large m, the numerical data (blue dots) falls of as to~ 4 , with the coefficient matching our analytic 
prediction given by F™^ (blue line) . Subtracting this leading order behaviour, we find that the remainder falls off as 

to~ 6 (purple squares and line) and then m -8 (brown diamonds and line), as expected. 



IV. DISCUSSION 



This paper extends the work of Paper I to the case of equatorial geodesic orbits in the Kerr spacetime. However, 
this only reflects a subset of the possible geodesic orbits in that case. In general, geodesies of the Kerr spacetime 
do not lie in the equatorial plane. While our calculation could be extended to cover the case of these more generic 
geodesies, we have chosen here to restrict ourselves to the case of equatorial motion and work with the significantly 
simpler expressions that ensue, leaving the more general case for future work. 

In our analysis, we have made use of scalar spherical harmonics which are not particularly well suited to the 
Kerr spacetime and gravitational perturbations. A more appropriate choice of basis functions may be the spheroidal 
harmonics; it may be more sensible to compute regularization parameters for a spheroidal harmonic basis. However, 
from a practical perspective, most existing numerical self-force calculations already make use of lower order versions of 
the expressions given here. The example clearly shows that these existing calculations gain significant improvements 
in accuracy from the use of spherical harmonic expansions. In this way, the end justifies the means: despite not 
being a natural choice, the use of spherical harmonics is a good practical choice. Nevertheless, an adaptation of our 
calculation to the spheroidal harmonic basis and to other gauges would make the results applicable in a much wider 
range of contexts. 

The Lorenz gauge metric perturbation equation on Kerr spacetime has yet to be been shown to be fully separable. 
It is likely that this would require the development of tensor spheroidal harmonics, which are as-yet unknown. In 
the absence of these, it has not been possible to test the validity of our electromagnetic and gravitational Z-modc 
regularization parameters. However, as in Paper I, deriving the expressions by independent methods gives us strong 
confidence in our results. Another check is to set a = 0, in which case the results agree with those of Paper I, which 
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we know to be correct. 

Note that the m-mode scheme is not affected by this issue as it is equally applicable to both Schwarzschild and Kerr 
spacetimes and is likewise equally as valid in the gravitational, electromagnetic and scalar cases. The only caveat 
is that the m-mode regularization parameters are only guaranteed to be correct for an effective source derived from 
a compatible approximation to the singular field. Since there is a large amount of flexibility in the effective source 
approach, if one chooses an incompatible singular field approximation, the regularization parameters here must be 
modified appropriately. 
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